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Abstract. In this paper we consider non-linear differential equations which 
are closely related to the generating functions of Frobenius-Euler polynomi- 
als. Prom our non-linear differential equations, we derive some new identities 
between the sums of products of Frobenius-Euler polynomials and Frobenius- 
Euler polynomials of higher order. 



1. Introduction 

Let u G C with u ^ 1. Then the Frobenius-Euler polynomials are defined by 
generating function as follows: 

(1) Fu{t, x) = ^T^e"* = V H,,{x I w)^, (see [1,2]). 

n=0 

In the special case, a; = 0, iJ„(0 | u) — Hn{u) are called the n-th Frobenius-Euler 
numbers (see [2]). 
By (1), we get 

(2) Hn{x\u)^^lj'\x''-^Hi{u) for n£Z+=NU{0}. 

Thus, by (1) and (2), we get the recurrence relation for Hn{u) as follows: 

f 1-u if n = 0, 

(3) Ho{u) = l, {Hiu) + 1) -Hn{u)=\ 

[0 II n > 0. 

with the usual convention about replacing H{u)^ by Hn{u) (see [2,10,12]). 
The Bernoulli and Euler polynomials can be defined by 

e — 1 ^ — ' n\ e"^ + 1 ^ — ' n\ 

n— n— 

In the special case, a; = 0, Bn{0) = Bn are the n-th Bernoulli numbers and £",1(0) — 
En are the n-th Euler numbers. 

The formula for a product of two Bernoulli polynomials are given by 

(4) 

Brn{x)Bnix)^}^\[ ]n+[ m -ir+' ^ ., B,n+n, 

^\\Zr J J I m + n — 2r [m + ny. 



where m + n > 2 and (™) = , „ = rnirn-i)-im-n+i) ^ r-^ 3- 

— \nJ n!{m— n)! n! \ L ' - 

1 
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From (1), wc note that Hn{x \ —1) = En{x). In [10], Nielson also obtained 
similar formulas for En{x)Em{x) and Em{x)Bn{x). 

In view point of (4), Carlitz have considered the following identities for the 
Frobenius-Euler polynomials as follows: 

(5) 

Hm{x I a)K^{x I /3) = | ap) ^^ ~ f^^^ ~ 

1 — ap 

+ '^^r^Y,{")Hr{a)Hr^+n-r{x \ a/3) + ^^^E {""^HsmHrn+n-six \ C 

r=0 '' s=0 

where a, /3 G C with a ^ 1, /3 7^ 1 and 7^ 1 (see [2]). 
In particular, if a 7^ 1 and aj3 = 1, then 



Hm{x I a)Hn(,x I = - (1 - a) ^ ()Hr{a) 



r=l 



r/ '^"'m + n — r + 1 



- (1 -a-i)y f"V^(""')-^"^""-'^'^'''^ 



,,s/ m + n — s + 1 

For r G N, the n-th Frobenius-Euler polynomials of order r are defined by generating 
function as follows: 

x) = F^{t, x) X F^{t, x) X ■ • • X a;) 



(6) 



r— times 

(1^) X (^^) X . . • X f e-* 

Ve* — u/ Ve* — u/ Ve* — M/ 



= H^^\x I u)— for u e C with u 7^ 1. 

In the special are called the n-th Frobenius-Euler 

numbers of order r (see [1-15]). 

In this paper we derive non-linear differential equations from (1) and wc study 
the solutions of non-linear differential equations. Finally, we give some new and 
interesting identities and formulae for the Frobenius-Euler polynomials of higher 
order by using our non-linear differential equations. 

2. Computation of sums of the products of Frobenius-Euler numbers 

and polynomials 



In this section we assume that 

F = F{t) = -fzr:^, and F^{t, x) = Fx ■•■ x F e^* for N eN. 

N— times 

Thus, by (7), we get 

at (e* — uy &■ — u (e* — uy 
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By (8), we get 

F^^\t,x) = F(i)(i)e*^ = -F{t,x)+uF^{t,x), 
and F^i) +F = uF'^. 
Let us consider the derivative of (8) with respect to t as follows: 

(10) 2uFF' =f" +f' . 
Thus, by (10) and (8), wc get 

(11) 2lu'^F^ -2uF^ =f" +f'. 
From (11), we note that 

(12) 2\u^F^ = f (2) + 3^ + 2F, where = 
Thus, by the derivative of (12) with respect to t, we get 

(13) 2!w23f2f' =i^(3)^3j.(2)^2FW, and F« = wF^ - F. 
By (13), we see that 

(14) 3\u^F^F = F^^^ + 6F(2) + IIF^^) + 6F. 
Thus, from (14), we have 

3!u*F*(t,a;) = F^^\t,x) + 6F^'^\t,x) + nF^^\t,x) + 6F{t,x). 
Continuing this process, we set 

N-l 

(15) {N - ly.u^-^F^ = ^ afe(iV)F('=), 

fe=0 

where F^*) = and TV e N. 

Now we try to find the coefficient ak{N) in (15). From the derivative of (15) 
with respect to t, we have 

Af-l N 

(16) iV!w^-iF^-iF(i) = ^ afc(iV)F('=+i) = ^ afc_i(iV)FW. 

k=0 fc=l 

By (8), we easily get 

(17) 7V!w^-1f^-1f(i) = 7V!n^-iF^-i(nF2 _ ^ An«^F^+i - N\u''-^F''. 
Prom (16) and (17), we can derive the following equation (18): 

N 

mu^F^+^ = N{N - l)\u^-^F^ + ak-i{N)F^''^ 

(18) ''^^ 

^ ' N-l N 

= TV ^ a,(JV)F« + ^ a,_i(7V)F«. 

fc=0 /s=l 

In (15), replacing N hy N + 1, we have 

TV 

(19) N\u^F^+^ =Yak{N + l)F^''\ 

fe=0 
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By (18) and (19), we get 

N 

(20) 

^ ' N-l N 

fe=0 fc=l 

By comparing coefficients on the both sides of (20), we obtain the following equa- 
tions: 

(21) Nao{N) = ao{N+l), ajv(iV + 1) = aiv-i(iV). 
For 1 < < n — 1, we have 

(22) ak{N+l) = Nak{N) + ak-i{N), 

where ak{N) = iov k > N or k < 0. From (21), we note that 

(23) ao{N + 1) = Nao{N) = N{N - l)ao{N - 1) = • • • = N{N -!)■■■ 2ao(2). 

By (8) and (15), we get 

1 

(24) F + f' =uF'^ ^Y,^k (P)F^^^ = ao {2)F + ai {2)F^^'> . 

k=0 

By comparing coefficients on the both sides of (24), we get 

(25) ao(2) = 1, and ai(2) = 1. 

From (23) and (25), we have ao{N) = [N — 1)!. By the second term of (21), we see 
that 

(26) aN{N + 1) = aN-i{N) = ajv-2(iV - 1) = • • • = ai(2) = 1. 

Finally, we derive the value of ak{N) in (15) from (22). 

Let us consider the following two variable function with variables s, t: 

(27) 9{t,s)=Y^ Yl «feW]v!*''' where |f|<l. 

iV>10<fc<JV-l 

By (22) and (27), we get 

N>lO<k<N-l 

fN j-N 

(28) =E E Nak^,{N+l)-s'^+J2 ^^^W^n 

Ar>10<fe<Af-l ■ Ar>10<fe<W-l 



Y E Nak+i{N)—s' + g{t,s). 

N>10<k<N-l 
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It is not difficult to show that 
(29) 

iV>10<fe<JV-l ■ JV>10<fe<JV-l 



{N-iy. 



4E E E -(^)(^ 

JV>ll<fe<JV ^ iV>l \0</c<JV ^ 

JV>1 \0<fe<Af ^ / \N>10<k<N ^ ' t 



= ;(.(M)- — 
From (28) and (29), we can derive the following equation: 

(30) E E «fe+i(^+l)-^ = ;(5(i>«)-Y— ^J+5(i,5). 



iV>l 0<fe<JV-l 

(31) 



'' k 



The left hand side of (13) = E E "'=+i(^)7;v — Tv* 

N>21<k<N-2 ^ '' 



JV>2 1<fe<Ar-l ^ ' \N>21<k<N-l 



s\^V ^ '(Af-1)! '(N-l)\ 

\N>2 0<k<N-l ^ ' ^ ' , 

\Ar>2 0<fe<JV-l ^ ^ / 

= ;(e E Miv)^.'-».(i)-^)=i(.(M)-^ 

By (30) and (31), we get 

(32) ,(M) + J(.'(M)-^)4(.(M)-^). 
Thus, by (32), we easily see that 

(33) = 5(t, s) + (<, s) + _ = <?(i, s) + <? (i, s) + -. 

By (33), we get 

(34) g(t,s) + ^-^<i{t,s) = --. 

s s 

To solve (34), we consider the solution of the following homogeneous differential 
equation: 

(35) o = 5(t,s) + ^5'(t,s). 
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Thus, by (35), we get 

(36) -g^t,s) = *-^g\t,s). 

By (33), we get 

g'{t,s) s 



^^^^ 9it,s) 1-f 

Prom (37), we have the following equation: 

(38) log5(i,s) = -slog(l-t) + C. 
By (38), we see that 

(39) 5(t,s) = e-"'°s(i-t)^ where A = e^. 
By using the variant of constant, we set 

(40) A = X{t, s). 
From (39) and (40), wc note that 

9it,s) = = A'(M)e---(-*) + '-Mf^'-^ 

(41) ,^ , 
= A'(t,s)e-^'°s(i-*) + ^s, 

where X'{t,s) ^ 

By multiply on both sides in (41), wc got 

(42) t_lg\t,s)+g{t,s) = A'^e-'^'^^^i-*). 
Prom (34) and (42), we get 

(43) _l = A'^e-«'°s(i-t). 

s s 

Thus, by (43), we get 

(44) a' =A'(t,s) = (l-i)'*-^ 
If we take indefinite integral on both sides of (44), we get 

(45) A = y Xdt = y (1 - ty-'^dt = -^(1 - ty + Ci, 

where Ci is constant. 

By (39) and (45), we easily see that 

(46) g{t, s) = e--^^°8(i-*) (--(l^ty+ Ci) . 

s 

Let us take t = in (46). Then, by (27) and (46), we get 

(47) o=--+Ci , Ci = -. 

s s 



-s 
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Thus, by (46) and (47), wc have 

g{t,s) = e-^>°g(i-') (---{1- ty) = -(1 - t)-^{i - (1 - ty) 

^ U - ^ ^ i / _^iog(i-t) _ 

s '' 
From (48) and Taylor expansion, we can derive the following equation (49): 

p(M) = iE^(-'»^(i-")" = SV St 

n>l n>l \ i=l / 



n>l \h = l i„=l 



Thus, by (49), we get 

^(*'«) = E(fcTT)! E ( E IIv^ 

fc>0 ^ ^ ' N>k+1 \h + -+lk+i=N 1^ '=+1/ 

Prom (27) and (50), we can derive the following equation (51): 

Therefore, by (15) and (51), we obtain the following theorem. 

Theorem 1. For u € C with u ^ 1, and N £ N, let us consider the following 
non-linear differential equation with respect to t: 

where F^^\t) = and F^{t) = F{t) x ■•■ x F{t). Then F{t) = ^ is a 

N— times 

solution of (52). 

Let us define F^''\t,x) = F^''\t)e*^. Then we obtain the following corollary. 

Corollary 2. For N eN, we set 

/V ^ 1 1 

(53) ^»,M) = ^i:^ E 

fe=0 iiH \-lk+i=N 

Then is a solution of (53). 
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From (1) and (6), we note that 



(54) 



e — w ni 

n=0 

and 

1 — UN / 1 — U\ / 1 — U 



(^)x(^)x-x(i^) = E-r(.4. 



n=0 



iV— times 



where Hn (u) are called the n-th Probenius-Euler numbers of order N. 
By (7) and (54), we get 



N— times 



1 /I — U\ / 1 - It \ 

~ il-u)^ Ve* - u) ^ Ve* - u) ^'"^ \e* -u) 

^ ^ ^ 

A''— times 

(55) 



^ ^ (=0 

and 

-(')=(H)(T^)-T^f:-'<4 



Prom (55), we note that 



(56) F''Ht) = '^=±m..iu)'^. 

1=0 



Therefore, by (52), (55) and (56), we obtain the following theorem. 
Theorem 3. For iV e N, n e Z+, we have 



fe=o ^ ^ ;i+-+ife+i=jv ^ ^ 
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From (55), wc can derive the following equation: 



N— times 



(57) 



oo 



V- / V- Hi,{u)Hi,iu)---Hi^{u)nl \t" 

n=0 \-In=0 i 1 jv j 



Therefore, by (57), we obtain the following corollary. 
Corollary 4. For N gN, n € Z+, we have 



E L ^J\HiAu)HiM---Hi,{u) 

JiH \-lN=n ^ ' ' ' 



By (53), we obtain the following corollary. 
Corollary 5. For G N, n G Z+, we /iaz;e 

, AT— 1 JV— 1 ^ 1 " / \ 

Prom (6), we note that 

/I — it\ /I — u\ / 1 — u 



n=0 



E~):,-(^)x(^)x-x(l^) 

(58) 



A^— times 



\n=0 / \)n=0 

n=0 \ 1=0 ' 



E E 

By comparing coefficients on both sides of (58), we get 
(59) H^^\x\u)=±(^'iy-^Hl''\u). 
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By the definition of notation, we get 



Coo J.! \ / °° m \ 

1=0 ■ / \ m=0 ■ / 



n=0 \ 1=0 



Prom (6), we note that 
(60) 

£»'S=(H)--(H) 



n=0 



N— times 

oo 



\/i=0 / \;,v=o / m=0 



n— \ iiH \-lN-\-m—n / 



By comparing coefficients on both sides of (58), we get 
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